function T=q1_T(x,y,D,t)
%计算每个镜子不同时间的单位转化矩阵

%已知参数
alpha_s=alp_s(D,t);%太阳高度角
gama_s=ga_s(D,t);%太阳方位角
z=ones(length(x),1)*4;%初始化所有镜子的高度
z0=8;%吸收塔高度

%初始化
T=cell(length(x),1);%用于存储每个镜子，镜面坐标系到地面坐标系转换单位矩阵



%thta_h为定日镜在镜场中的方位角 cos_thta为镜面中心主光线的入射余弦
for i=1:length(x)
    lambda(i)=atan((x(i)^2+y(i)^2)^0.5/(z0-z(i)));
    if y(i)>=0
        thta_h(i)=acos(x(i)/((x(i)^2+y(i)^2))^0.5);
    else
        thta_h(i)=atan(y(i)/x(i)+pi);%y<0的情况
    end
    cos_thta(i)=(2^0.5/2)*(sin(alpha_s)*cos(lambda(i))-cos(thta_h(i)-gama_s)*cos(alpha_s)*sin(lambda(i))+1)^0.5;
    
    %镜面中心法线的方位角
    A_H(i)=atan((sin(thta_h(i))*sin(lambda(i))-sin(gama_s)*cos(alpha_s))/(cos(thta_h(i))*sin(lambda(i))-cos(gama_s)*cos(alpha_s)));
    %镜面中心法线的俯仰角
    E_H(i)=acos((sin(alpha_s)+cos(lambda(i)))/(2*cos_thta(i)));
    
    %单个镜子的单位转换矩阵
    Temp=ones(3,3);
    Temp(1,:)=[-sin(E_H(i)) -sin(A_H(i))*cos(E_H(i)) cos(A_H(i))*cos(E_H(i))];
    Temp(2,:)=[cos(E_H(i)) -sin(A_H(i))*sin(E_H(i)) cos(A_H(i))*sin(E_H(i))];
    Temp(3,:)=[0 cos(A_H(i)) sin(A_H(i))];
    
    T{i,1}=Temp;
end

 

end

